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Abstract 

The qq spectrum is studied in a generalized constituent quark model con- 
strained in the study of the AW phenomenology and the baryon spectrum. 
An overall good fit to the available experimental data is obtained. A de- 
tailed analysis of all sectors from the light-pseudoscalar and vector mesons 
to bottomonium is performed paying special attention to the existence and 
nature of some non well-established states. These results should serve as a 
complementary tool in distinguishing conventional quark model mesons from 
glueballs, hybrids or multiquark states. 
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I. INTRODUCTION 



Meson spectroscopy is an extremely broad subject covering from the few hundred MeV 
masses of the light pseudoscalar mesons to the 10 GeV scale of the bb system. Such a wide 
energy region allows to address perturbative and nonperturbative phenomena of the under- 
lying theory. The continuously increasing huge amount of data and its apparent simplicity 
made mesons ideal systems to learn about the properties of QCD. As an exact solution of 
the theory seems not attainable at present and the big effort done in developing its lattice 
approximation is still not able to produce, without guidance, realistic results, phenomeno- 
logical models are the most important theoretical tool to study the meson properties. Last 
years are being extremely exciting due to the new data obtained at the B factories, reopen- 
ing the interest of classifying the meson experimental data as qq states according to SU (N) 
irreducible representations. The recently measured -D*j(2317) and X(3872) states seem to 
be hardly accommodated in a pure qq description. The D*j(2317) presents a mass much 
lower than the prediction of potential models, whereas the X(3872) is too light to be a 2P 
charmonium state and too heavy for a ID state. Such discrepancies could be related with 
the structure of some of the light scalar mesons still not well established theoretically. To 
be able to understand the nature of new resonances it is important that we have a template 
against which to compare observed states with theoretical predictions. 

Much has been learned during the last years about the structure and properties of QCD. 
The study of charmonium and bottomonium made clear that heavy-quark systems are prop- 
erly described by nonrelativistic potential models reflecting the dynamics expected from 
QCD [1]. The a priori complicated light-meson sector was quite surprisingly well reproduced 
on its bulk properties by means of a universal one-gluon exchange plus a linear confining 
potential [2]. However, the dynamics of the light-quark sector is expected to be dominated 
by the nonperturbative spontaneous breaking of chiral symmetry, basic property of the QCD 
Lagrangian not satisfied by the first constituent quark model approaches to the light-meson 
sector. Chiral symmetry breaking suggests to divide the quarks into two different sectors: 
light quarks (u, d, and s) where the chiral symmetry is spontaneously broken, and heavy 
quarks (c and b) where the symmetry is explicitly broken. The origin of the constituent 
quark mass can be traced back to the spontaneous breaking of chiral symmetry and conse- 
quently constituent quarks should interact through the exchange of Goldstone bosons [3], 
these exchanges being essential to obtain a correct description of the NN phenomenology 
and the light baryon spectrum. Therefore, for the light sector hadrons can be described as 
systems of confined constituent quarks (antiquarks) interacting through gluons and boson 
exchanges, whereas for the heavy sector hadrons are systems of confined current quarks 
interacting through gluon exchanges. Concerning the other basic property of the theory, 
confinement, little analytical progress has been made. Lattice calculations indicate that the 
linear confining potential at short distances is screened due to pair creation making it flat 
at large distances [4]. 

A preliminary study of the light-meson spectra and their strong decays has already been 
done with an interacting potential constrained on the study of the NN system and the 
nonstrange baryon spectrum [5]. Our purpose is to generalize this work obtaining a quark- 
quark interaction that allows for a complete study of the meson spectra, from the light to 
the heavy sector. To find new physics, it is important that we test the quark model against 
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known states to understand its strengths and weakness. To this end we shall begin in the 
next section discussing the theoretical ingredients of the constituent quark model. In Sec. 
Ill we perform a detail comparison of the predictions of our model with experiment from 
the light pseudoscalar and vector mesons to bottomonium. This will allow us to identify 
discrepancies between the quark model predictions and experiment that may signal physics 
beyond conventional hadron spectroscopy. We shall go over these puzzles to decide whether 
the discrepancy is most likely a problem with the model or the experiment, or whether 
it most likely signals some new physics. In particular the possible presence of four quark 
systems in the scalar meson sector will be addressed. Finally, in Sec. IV we summarize our 
most important findings. 



II. SU(3) CONSTITUENT QUARK MODEL 

Since the origin of the quark model hadrons have been considered to be built by con- 
stituent (massive) quarks. Nowadays it is widely recognized that the constituent quark 
mass appears because of the spontaneous breaking of the original SU(3)l <S> SU(3)r chiral 
symmetry at some momentum scale. The picture of the QCD vacuum as a dilute medium 
of instantons [6] explains nicely such a symmetry breaking, which is the most important 
nonperturbative phenomenon for hadron structure at low energies. Quarks interact with 
fermionic zero modes of the individual instantons in the medium and therefore the propaga- 
tor of a light quark gets modified and quarks acquire a momentum dependent mass which 
drops to zero for momenta higher than the inverse of the average instanton size p. The 
momentum dependent quark mass acts as a natural cutoff of the theory. In the domain of 
momenta k < 1/p, a simple Lagrangian invariant under the chiral transformation can be 
derived as [6] 

L = ^(i 7 ^ - MIP 6 )^ (1) 

being U 15 = exp(i7r a A a 75// 7r ). 7r a denotes the pseudoscalar fields (n,Ki,ri s ) with i=l,...,4, 
and M is the constituent quark mass. An expression of the constituent quark mass can be 
obtained from the theory, but it also can be parametrized as M = m q F(q 2 ) with 



f(q 2 ) 



A 2 



A 2 + q 2 



(2) 



where A determines the scale at which chiral symmetry is broken. Once a constituent quark 
mass is generated such particles have to interact through Goldstone modes. Whereas the 
Lagrangian ^(ry^c^— M)^ is not invariant under chiral rotations, that of Eq. (1) is invariant 
since the rotation of the quark fields can be compensated renaming the bosons fields. £/ 75 
can be expanded in terms of boson fields as, 

t/ 75 = 1 + ^7 5 A a 7T a - -^TtV* + ... (3) 

The first term generates the constituent quark mass and the second one gives rise to a one- 
boson exchange interaction between quarks. The main contribution of the third term comes 
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from the two-pion exchange which can be simulated by means of a scalar exchange potential. 
Inserting Eqs. (2) and (3) in Eq. (1), one obtains the simplest Lagrangian invariant under 
the chiral transformation SU (3) l <S> SU (3)r with a scale dependent constituent quark mass, 
containing SU (3) scalar and pseudoscalar potentials. The nonrelativistic reduction of this 
Lagrangian has been performed for the study of nuclear forces and will not be repeated 
here, although the interested reader can follow the particular steps in several theoretical 
works [7,8]. The different terms of the potential contain central and tensor or central and 
spin-orbit contributions that will be grouped for consistency. Therefore, the chiral part of 
the quark-quark interaction can be resumed as follows, 



= vfJnj) + vlfc) + v q s q °(f tJ ) , 



(4) 



where C stands for central, T for tensor, and SO for spin-orbit potentials. The central part 
presents four different contributions, 
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the angle 6p appears as a consequence of considering the physical rj instead the octet one. 
g c h = m q j f n , the A's are the SU(3) flavor Gell-Mann matrices, is the quark mass and m n , 
mx and m n are the masses of the SU(3) Goldstone bosons, taken to be their experimental 
values. m a is determined through the PCAC relation m 2 ~ m 2 + Am 2 ud [9]. Finally, Y(x) 
is the standard Yukawa function defined by Y(x) = e~ x fx. 

There are three different contributions to the tensor potential, 
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each term given by, 
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being Sij = 3 (<7j • rij)(cfj • fij) — &i • &j the quark tensor operator and H(x) = (1 + 3/x + 
3/x 2 )Y(x). 

Finally, the spin-orbit potential only presents a contribution coming form the scalar part 
of the interaction, 



Ait A 2 — m 2 2rriimj 



A 3 

G(m rij) - -%G{Knj) 



L-S (9) 



where G(x) = (1 + l/x)Y(x)/x. The chiral coupling constant is determined from the 
tiNN coupling constant through 

9ch _ (3\ 2 9lNN m u,d /-, n \ 

what assumes that flavor SU (3) is an exact symmetry only broken by the different mass of 
the strange quark. 

This interaction, arising as a consequence of the instanton induced chiral symmetry 
breaking, gives rise, among other effects, to vector-pseudoscalar meson mass splitting and it 
also generates flavor mixing for the 77 mesons. This is the very same effect obtained by other 
instanton induced approaches as those based on the 't Hooft Lagrangian [10,11]. In the 
heavy-quark sector chiral symmetry is explicitly broken and the Goldstone boson exchange 
interaction is not active in such a way that we cannot reproduce the hyperfme splittings 
for heavy mesons. This is a common feature of other instanton induced interactions [12] 
unless a phenomelogical parametrization is done [13]. Beyond the chiral symmetry breaking 
scale one expects the dynamics being governed by QCD perturbative effects. They mimic 
the gluon fluctuations around the instanton vacuum and are taken into account through the 
one-gluon-exchange (OGE) potential. Such a potential nicely describes the heavy-meson 
phenomenology. Following de Rujula et al. [14] the OGE is a standard color Fermi-Breit 
interaction given by the Lagrangian, 

L = iV&ra s ^G>*\ c ip (11) 

where A c are the SU(3) color matrices, G^ is the gluon field and a s is the quark-gluon 
coupling constant. The nonrelativistic reduction of the OGE diagram in QCD for point-like 
quarks presents a contact term that, when not treated perturbatively, leads to collapse [15]. 
This is why one maintains the structure of the OGE, but the S function is regularized in 
a suitable way. This regularization is justified based on the finite size of the constituent 
quarks and should be therefore flavor dependent [16]. As a consequence, the central part of 
the OGE reads, 

V OGE {r i3 ) = -a. X c i ■ X c j j- - —a t • a 3 j (12) 

where r (fi) = f //i, scaling with the reduced mass as expected for a coulombic system. Let 
us note that the nonrelativistic reduction of the one-gluon exchange diagram in QCD yields 
several spin-independent contributions [14] that have not been considered. As the OGE is 
an effective interaction we have included only the relevant different structures obtained from 
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the nonrelativistic reduction and neglected the other terms that are supposed to be included 
in the fitted parameters [2]. 

The noncentral terms of the OGE present a similar problem. For point-like quarks they 
contain an 1/r 3 term that in spite of its strength has been usually treated perturbatively. 
Once again the finite size of the constituent quarks allows for a regularization and therefore 
for an exact treatment of these contributions, obtaining tensor and spin-orbit potentials of 
the form, 



VoGE^fij) — 



a. 



16 mirrij 1 
16 m?m? * 3 



1 
1 

ij 



,-rij/r g (n) 



~ r ij/r g (p) 



+ 



1 + 



3r» 



+ 



X 



{{mi + rrij) 2 + 2m i m j ) (S + ■ L) + (m) - rrif) (£_ ■ L) 



(13) 



where S± — Si ± Sj, and r s (/x) = r g /pL presents a similar behavior to the scaling of the 
central term. The wide energy covered to describe the light, strange and heavy mesons 
requires an effective scale-dependent strong coupling constant [17] that cannot be obtained 
from the usual one-loop expression of the running coupling constant because it diverges 
when Q — > Aq C d- The freezing of the strong coupling constant at low energies studied in 
several theoretical approaches [18,19] has been used in different phenomenological models 
[20]. The momentum-dependent quark-gluon coupling constant is frozen for each flavor 
sector. For this purpose one has to determine the typical momentum scale of each flavor 
sector that, as explained in Ref. [21], can be assimilated to the reduced mass of the system. 
As a consequence, we use an effective scale-dependent strong coupling constant given by 



a s (n) 



a 



(14) 



where fi is the reduced mass of the qq system and ao, /xo an d A are determined as explained 
in Sec. III. This equation gives rise to a s ~ 0.54 for the light-quark sector, a value consistent 
with the one used in the study of the nonstrange hadron phenomenology [5,22], and it also 
has an appropriate high Q 2 behavior, a s ~ 0.127 at the Z mass [23]. In Fig. 1 we compare 
our parametrization to the experimental data [24,25]. We also show for comparison the 
parametrization obtained in Ref. [18] from an analytical model of QCD. 

When Goldstone-boson exchanges are considered together with the OGE, the possibility 
of double counting emerges. This question connects directly with the nature of the pion, 
studied for a long time concluding its dual character as qq pair and Goldstone boson [26,27]. 
In first approximation, it should be reasonable to construct a theory in which chiral sym- 
metry is retained in the Goldstone mode but the internal structure of the pion is neglected. 
This would be essentially a long- wave length approximation [27]. This may be the reason 
why the OPE generates contributions that are not obtained from the OGE, while the p and 
uj meson exchanges give rise to contributions already generated by the OGE [28]. Explicit 
studies on the literature about the double counting problem concluded that while the pion 
can be safely exchanged together with the gluon, the vector and axial mesons cannot [29]. 
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Finally, any model imitating QCD should incorporate another nonperturbative effect, 
confinement, that takes into account that the only observed hadrons are color singlets. It 
remains an unsolved problem to derive confinement from QCD in an analytic manner. The 
only indication we have on the nature of confinement is through lattice studies, showing 
that qq systems are well reproduced at short distances by a linear potential. Such potential 
can be physically interpreted in a picture in which the quark and the antiquark are linked 
with a one-dimensional color flux tube. The spontaneous creation of light-quark pairs may 
give rise to a breakup of the color flux tube. It has been proposed that this translates into 
a screened potential [4], in such a way that the potential does not rise continuously but it 
saturates at some interquark distance. Although string breaking has not been definitively 
confirmed through lattice calculations [30], a quite rapid crossover from a linear rising to 
a flat potential is well established in SU(2) Yang-Mills theories [31]. A screened potential 
simulating these results can be written as, 



where A is a global constant fixing the origin of energies. At short distances this potential 
presents a linear behavior with an effective confinement strength a = a c fi c \ c i ■ \ c j, while 
it becomes constant at large distances. Such screened confining potentials provide with an 
explanation to the missing state problem in the baryon spectra [32] and also to the deviation 
of the Regge trajectories from the linear behavior for higher angular momentum states [33]. 

One important question which has not been properly answered is the covariance property 
of confinement. While the spin-orbit splittings in heavy-quark systems suggest a scalar 
confining potential [34], Ref. [35] showed that the Dirac structure of confinement is of vector 
nature in the heavy-quark limit of QCD. On the other hand, a significant mixture of scalar 
and vector confinement has been used to explain the decay widths of P-wave D mesons [36]. 
Nonetheless, analytic techniques [37] and numerical studies using lattice QCD [38] have 
shown that the confining forces are spin independent apart from the inevitable spin-orbit 
pseudoforce due to the Thomas precession [39]. Therefore, we will consider a confinement 
spin-orbit contribution as an arbitrary combination of scalar and vector terms, 



where a s would control the ratio between them. 

Once perturbative (one-gluon exchange) and nonperturbative (confinement and chiral 
symmetry breaking) aspects of QCD have been considered, one ends up with a quark-quark 
interaction of the form (from now on we will refer to a light quark, u or d, as n, s will be 
used for the strange quark and Q for the heavy quarks c and b): 
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The corresponding qq potential is obtained from the qq one as detailed in Ref. [40]. In 
the case of Vx(^j), where G parity is not well defined, the transformation is given by 
A"- A| — > A" - (A2) 71 , which recovers the standard change of sign in the case of the pseudoscalar 
exchange between two nonstrange quarks. 

Apart from models incorporating the one-gluon exchange potential for the short-range 
part of the interaction, one finds in the literature other attempts to study hadron phe- 
nomenology based on instanton induced forces. In Ref. [41] a flavor antisymmetric quark- 
quark instanton induced interaction was derived. It was used in a nonrelativistic framework 
for the study of the two-baryon system and the baryon spectrum, but it was never applied 
to study the meson spectrum. Such a model does not contain Goldstone boson exchanges, 
which are essential to make contact with the one and two baryon systems. The short-ranged 
instanton induced force was supplemented by a baryonic meson exchange potential to give a 
quantitative explanation of experimental data. Moreover, in Ref. [42] it was demonstrated 
that this instanton induced force reproduces the baryon spectrum as well as the one-gluon 
exchange. The other extensive work using instanton induced interactions is that of the 
group of Bonn [43,44]. The interaction is derived in Ref. [43] obtaining a sum of contact 
terms that are regularized (as we have done for the OGE potential) and whose flavor mixing 
matrix elements are fitted to experimental data. This interaction was supplemented by a 
phenomenological confining potential and applied in a Bethe-Salpeter framework to study 
the light meson and baryon spectra. This instanton induced force is only valid for light 
quarks in such a way that the extension of the model to study heavy flavors is done in a 
completely phenomenological way [13]. Therefore, these instanton induced models do not 
allow for the moment for a coherent study of the light and heavy flavors neither for a simul- 
taneous description of the baryon-baryon phenomenology. These two approaches, the one 
based on the 't Hooft interaction and the one followed in the present work, clearly shares 
the most important features of the quark-quark interaction. The reason why we adopt the 
first scheme lies in the fact that it includes explicitly Goldstone boson exchanges between 
quarks which are essential to make contact with the one and two-baryon phenomenology. 

III. RESULTS AND DISCUSSION 

Let us first discuss how the parameters of the model are fixed. Most of them, as for 
example those of the Goldstone boson fields, are taken from calculations on two-baryon 
systems. Once the Goldstone boson exchange part of the interaction has been determined, 
the one-gluon exchange controls the hyperfine splittings. It is worth to notice the relevance 
of the OGE contribution for the description of the meson spectra, as has already been 
emphasized in Ref. [5]. The conclusion of this work, that the p — 71 mass difference cannot 
be obtained from the Goldstone-boson exchanges alone is fully maintained. For the sake 
of completeness let us mentioned that in the present model this mass difference would be 
around 38 MeV if the OGE is not used. As it is also illustrated in Ref. [5] the dependence 
of the p — ii mass difference on the Goldstone-boson cut-off masses is small, in such a 
way that the spectra come never determined by the election of the cut-off masses. We 
have fixed the OGE parameters by a global fit to the hyperfine splittings well established 
by the Particle Data Group (PDG) [45] from the light to the heavy-quark sector. The 
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strength of confinement determines the energy difference between any J PC meson ground 
state and its radial excitations. We have fitted a c and fi c in order to reproduce two well 
measured energy differences: the p meson and its first radial excitation and the J/ip and 
the ip(2S), obtaining values close to those inferred in Ref. [46] through the analysis of the 
screening of QCD suggested by unquenched lattice calculations in the heavy-quark sector. 
Finally, it only remains to fix the relative strength of scalar and vector confinement. For this 
purpose one has to look for some states where the spin-orbit contribution, being important, 
can be easily isolated from other effects as for example flavor mixing. This is the case 
of the isovector mesons ai(1260) and a 2 (1320). In Fig. 2 we have plotted their masses 
as a function of a s , the relative strength of scalar and vector confinement. As can be 
seen, using a strict scalar confining potential, a s — 1, one would obtain 1343 MeV and 
1201 MeV, respectively, in complete disagreement with the ordering and magnitude of the 
experimental data. Introducing a mixture of vector confinement, a s =0.777, the experimental 
order is recovered being now the masses 1205 MeV and 1327 MeV, respectively, both within 
the experimental error bars. This value also allows to obtain a good description of the 
experimental data in the cc and bb systems. The parameters of the model are resumed in 
Table I. 

Acceptable results for the meson spectra have been provided by relativistic as well as 
nonrelativistic approaches [48]. In both cases a QCD-inspired interaction is used, the dif- 
ference being in quarks masses and kinematics. In fact, several works in the literature [49] 
have studied the connections existing between relativistic, semirelativistic, and nonrelativis- 
tic potential models of quarkonium using an interaction composed of an attractive Coulomb 
potential and a confining power-law term. The spectra of these very different models be- 
come nearly similar provided specific relations exist between the dimensionless parameters 
peculiar to each model. 

As a consequence we will solve the Schrodinger equation for the relative motion of the 
qq pair with the interacting potential of Eq. (17). There is a particularly simple and 
efficient method for integrating this type of second order differential equations, the commonly 
called Numerov method [50]. The noncentral potentials (tensor and spin-orbit) give their 
most important contribution for the diagonal terms. For example for the isovector states 
J PC = 1 (L = or 2, S — 1, J = 1) we would have the following matrix elements (in 
MeV): 

/ < 3 S'xI^Tl 3 ^! >= 772 < 3 S' 1 | J f/| 3 Z? 1 >= 22 \ 

V < 3 .Dil^ri 3 ^! >= 22 < 3 D l \H\ 3 D l >= 1518/ ' ^ } 

and for the isovector J PC = 2++ (L = 1, 3, S = 1, J = 2) 

(19) 



< 3 P 2 \H\ 3 P 2 >= 1327 < 3 P 2 \H\ 3 F 2 >= 5 



2 \n\ r 2 ^— j-ozi ^ r 1\ n \ r 2 

< 3 F 2 \H\ 3 P 2 >= 5 < 3 F 2 \H\ 3 F 2 >= 1797 

where the matrix elements are calculated with the wave functions solution of the single- 
channel Schrodinger equation for the total hamiltonian, including central and noncentral 
terms. The perturbative effect of the nondiagonal contributions of the noncentral potentials 
is observed. The importance of the diagonal contributions can be easily inferred from the 
energies obtained when the different interactions are connected. In table II we give the 
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energies obtained when solving the Schrodinger equation for two different partial waves: one 
of the cases discussed above, the isovector 3 Di, and the isoscalar 3 P (rin). One can see that 
the corrections induced by the diagonal contributions of the noncentral potentials are much 
important than the nondiagonal ones. This effect becomes clearly nonperturbative for the 
isoscalar 3 P state, that on the other hand does not admit coupling to partial waves with 
different orbital angular momentum. Therefore the noncentral potentials will be treated 
exactly in their diagonal contributions and by diagonalizing the corresponding matrix when 
nondiagonal contributions are present. The same reasoning as before applies in the case of 
the spin-orbit contribution proportional to S- as will be discussed in Sect. 111H. 

In Tables III to XII we compare the masses obtained within the present model to ex- 
perimental data [45]. Tables III, IV, and V do not include the scalar mesons which will be 
discussed separately in Sec. Ill L and are shown in Table XII. In all cases we show the name 
and mass of the experimental state and the prediction of our model together with the J PC 
quantum numbers, or J p if C parity is not well defined (in these cases we have explicitly 
indicated the spin). If there is an experimental meson without assignment of quantum num- 
bers, those indicated by question marks, and we find a corresponding state in our model, 
we indicate between square brackets its quantum numbers. In the tables we also include in 
parenthesis the radial excitation and the orbital angular momentum corresponding to the 
state under consideration. 

A final comment about the flavor mixing angles is in order. Our interacting hamiltonian 
needs as input the pseudoscalar octet-singlet mixing angle (see Eq. (6)). For this angle, 
values in the range of —10° to —23° have been obtained depending on the analysis performed 
[45,51-56]. We have taken an intermediate value of —15°. The model provides a theoretical 
mixing nn <-> ss (uj <-> 0, i] <-> 77', ...) as a consequence of the Vk potential. The mixing 
angle obtained is 9 V = 34.7° for the u <-> case, and 9 P = —21.7° for the 77 <-> rj . We want 
to emphasize the independence of our result on the input value for dp. For example if we 
had used a theoretical input of 9p = —23°, our predicted pseudoscalar mixing angle would 
have been dp = —21.9°, almost the same as before. Our result is compatible with many 
others reported in the literature [52,54,56] and close to the 6 l p n value given on the PDG 
[45] (6 l p n = —23°), as expected because we use a mixing formula based on linear and not 
quadratic masses. The mixing angle obtained for most meson nonets is approximately ideal, 
exceptions are the nonet of pseudoscalar and scalar mesons. The mixing angle, calculated 
for the ground state of each J PC nonet, is assumed for the radial excitations with the same 
quantum numbers. 

There are a number of states whose quantum numbers are not clearly determined or do 
not seem to present a clear correspondence with a qq state. Let us analyze in detail the 
results. 

A. / = J PC = 0-+ states 

There is an overall good agreement between the constituent quark model results and the 
experimentally observed mesons. Recently there has been a modification in the experimental 
situation. In the 2002 PDG [71] there were two resonances in the 1.5 GeV energy region, 
77(1295) and 77(1440). However, the 2004 PDG reports three states, 77(1295), 7/(1405) and 
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77(1475), all of them interpreted as pseudoscalar h states. This modification was mainly 
due to a recent experiment by the E852 Collaboration [72]. They found evidence for three 
pseudoscalar resonances: 7/(1295), 77(1416) and r/(1485) in the analysis of the reaction n~p — > 
K + K~ir°n. The first one decays exclusively into a 7r°, the second into a 7r° and K*K and 
the last one exclusively into K*K. We only obtain two CT + states in this energy region, 
one at 1290 MeV and the other at 1563 MeV, both being predominantly nn. Based on our 
results one would identify the first state with the 7/(1295). However, our model predicts 
only one state corresponding either to 7/(1405) or 7/(1475). Although the energy seems 
to favor an assignment of our second state with the 7/(1475), its dominant nn component 
makes difficult to explain its experimental decay modes (K*K). However, the expectation 
of the 7/(1475) being a dominantly ss 2 1 Sq (sometimes the discussion is more clear using 
the spectroscopic notation n 2S+1 Lj) state [73] as suggested by the observation of a large 
K*K decay mode is highly dependent on the strong KK final state interaction [74]. The 
existence of the third resonance would therefore imply the presence of additional states 
beyond the two obtained in the quark model. However the experimental situation is not 
definitively settled. There are some speculations that one of these states could be somehow 
related to the /i(1420) [75]. To disentangle the flavor content of the 7/(1405) and 7/(1475) 
could be a very important experimental contribution feasible at CLEO, what would help to 
discriminate among different theoretical models. 

To illustrate the uncertain situation with these resonances let us finally mention that 
the analysis of the KKn and nmr channels in 77 collision performed in Ref. [76] observed 
the 7/(1475) in KKn, but not the 7/(1405) in nmr. Since gluonium production is presumably 
suppressed in 77 collisions this result suggests that the 7/(1405) may have a large glueball 
component [77]. This interpretation, however, is not favored by lattice calculations, which 
predict the 0~ + state above 2.5 GeV. 

B. / = J PC = 1++ states 

The historically confused experimental status of light axial vectors has improved a lot 
with high statistic central production experiments on nim, K Kn and 4ir by WA102 [78] 
and K Kit by E690 [79]. In these experiments very clearly /i(1285) and /i(1420) states have 
been observed, but there is no evidence of /i(1510). Although the /i(1420) and the /i(1510) 
are well separated in mass and resolved in different experiments, there are no experiment 
or production reactions in which both states have been detected. All these observations 
express skepticism regarding the existence of the /1 (1510) [80]. In view of their masses the 
obvious assumption is that the /i(1285) is the 1 3 Pi nn state and the /i(1420) its ss partner. 
In Ref. [80] a significant singlet-octet mixing for the nonet of the /i(1285) and /i(1420) has 
been obtained from several independent analysis, 9 ~ 50°. We obtain a value of 46.3° in 
complete agreement with the former study. One should remember that the mixing in our 
model is fully determined by the structure of the potential. 
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C. / = J PC = 2~+ states 

Regarding the 1 = 2~ + the situation seems to be clear. Experimentally there are 
two states: the ?7 2 (1645) with a mass of 1617±5 MeV and decaying into a 2 (1320)7r, and 
the 7/ 2 (1870) with a mass of 1842±8 and decaying into a 2 (1320)7r and Art. We obtain three 
theoretical states in this energy region: the first one with a mass of 1600 MeV and pure 
light-quark content, a second with a mass of 1853 MeV and pure ss content and a third 
state with a mass of 1863 MeV and also pure light-quark content. The assignment of our 
first state to the ?7 2 (1645) seems clear. With respect to the ?7 2 (1870) its reported decay 
modes are not accessible in the 3 P model having a pure ss content. This enforces to assign 
the ?7 2 (1870) to the second excited state with light-quark content. Thus, it may exist an 
unobserved resonance close to the ?7 2 (1870) being a pure ss state. 

D. / = J PC = 1— and J PC = 3~ states 

In the PDG there are three isoscalar 1 states in the 1.5 GeV energy region: u;(1420), 
o»(1650) and 0(1680). The first and third state seem to be fairly well established, however 
in the last two years there have been several modifications on the mass of the o>(1650). 
The 2002 PDG gave a mass of 1649 ± 24 MeV, the 2003 electronic version reported a mass 
between 1600 - 1800 MeV, while the 2004 PDG quotes a value of 1670±30 MeV. This mod- 
ification is based on the analysis of the reaction e + e~ — > [81] where two uj states are 
reported, the first one with a mass of 1490 ± 75 MeV decaying into 7r + 7r~7r° and the second 
with a mass of 1790 ± 50 and decaying, with approximately the same probability, into 37r 
and ujTTTT . While the first state clearly corresponds to the cu(1420), the second state was 
included by the electronic 2003 PDG as an experimental result for the u;(1650), increasing 
drastically the error bars. In the 2004 PDG this experiment has not been considered within 
the statistical fit, therefore reducing again the error bars. A similar situation may happen 
with the 0(1680). The FOCUS Collaboration at Fermilab [82] has reported a high statis- 
tic study of the diflractive photoproduction of K + K~ confirming the existence of a clear 
enhancement of the cross section corresponding with a fitted mass of 1753 ± 4 MeV. 

Theoretically we find three 1 states in this energy region: two almost pure light-quark 
states with masses of 1444 MeV and 1784 MeV, and an almost pure ss state with a mass 
of 1726 MeV. Based on the aforementioned discussion, the first two states may correspond 
to the a;(1420) and c<j(1650), whereas the third one would correspond to the 0(1680). Our 
analysis also assigns the cj(1650) to a 3 3 Si state instead of a l 3 Di wave [73], that would 
have a mass of 1475 MeV with a pure light-quark content, in complete disagreement with 
the experimental data. 

For the 3 case we obtain two almost degenerate states with completely different flavor 
content. Experimentally there is a clear evidence that the 03(1850) has a dominant ss flavor 
content, and therefore it would correspond to our ss state at 1875 MeV, predicting the 
existence of an unobserved nn 3 state in the same energy region. 
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E. / = J PC = 2++ states 

Experimentally there is a proliferation of isoscalar 2 ++ states in an energy region that 
has been suggested as coexisting with 2 ++ glueballs. Theoretically the results of our model 
may confirm that the / 2 (1270) and /2(1525) are the nn and ss members of the 1 3 P 2 <?<? flavor 
nonet. Although nn <-> ss mixing is present due to the kaon exchange, the nn content of 
the /2(1525) is smaller than 0.1% in agreement with the experimental /2(1525)77 coupling, 
which limits the nn content to a few percent [74]. Our model also finds that the / 2 (1565), 
observed in pp annihilation at rest, and the / 2 (1640), decaying into ujuj and 4tt, could be 
the same state, the nn member of the 2 3 P 2 qq flavor nonet as is suggested in the PDG. 
The / 2 (1950) would be its ss partner in agreement with the experimentally observed decays. 
Finally the / 2 (1810) corresponds to the nn member of the 1 3 P 2 nonet and the / 2 (1910) to 
the nn member of the 3 3 P 2 nonet, as expected from the decay patterns. 

In this sector there is a meson, the / 2 (1430), that has no equivalent in our qq scheme. This 
state is not confirmed in the PDG and even recent measurements have suggested a different 
assignment of quantum numbers, being a ++ state [83], what could make it compatible with 
the lightest scalar glueball [84]. It was seen together with another resonance, the / 2 (1480), in 
an experiment where the / 2 (1270) was not detected [85]. A full understanding of the nature 
of the / 2 (1430) will probably also require an explanation for the absence of the / 2 (1270) in 
the experimental data of Ref . [85] . 

F. / = J PC = 1+- states 

While the nn 1 3 Pi quark model prediction, 1257 MeV, differs slightly from its corre- 
sponding experimental state, /i(1285), the nn l l P\ is a little bit far from the expected 
corresponding experimental state, /ii(1170). Annihilation contributions could improve the 
agreement with data, having in mind that they are expected to be negative and larger in 
the S = light-quark sector while almost negligible for S — 1 states [2]. There are some 
estimations of the contribution of annihilation in the light-quark sector, but its dependence 
on fitted parameters prevents from making any definitive conclusion. The /i!(1380) it is a 
convincing candidate for the ss partner of the 1 l P\ /ii(1170). Although a first glance to the 
quark model result seems to indicate that the mass is too high, the most recent experimental 
measurement for this state reports a mass of 1440±60 [86]. The last measured l + ~ state, the 
/ii(1595), will correspond to the nn 2 l P\ state, its ss partner being an unobserved meson 
around 1973 MeV as also noticed from the analysis of strange decays in Ref. [74]. 

G. / = 1 mesons 

The results shown in Table IV present an almost perfect parallelism between the theoreti- 
cal predictions and the experimentally observed states. The only not well defined correspon- 
dence comes from the p(1700), which is compatible with an S or D wave. Our calculation 
indicates a mass for this state a little bit higher than the one reported in the PDG and 
compatible with the observations of Ref. [87]. 
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H. 1 = 1/2 strange mesons 



From a theoretical perspective strange mesons exhibit explicit flavor and therefore they 
do not present the additional complications of annihilation mixing what makes the isoscalar 
mesons a priori much more complicated. 

There is a good correspondence between the 0~ and 1~ quark model results and the 
experimental mesons. In particular, the P(1630), whose quantum numbers have not been 
yet determined, would correspond to the 2 3 Si ns state. The only problem in this sector 
appears for the if* (1410), in clear disagreement with the predictions of the quark model. 
The if* (1410) could be the radial excitation of the if* (892), because it is the first 1" vector 
resonance observed. However, the first ns radial excitation is predicted at 1620 MeV. This 
is reasonable taking into account that in the 5=1 light-quark sector the radial excitation 
has a mass around 700 MeV above the ground state, what confirms that the first radial 
excitation of the K * (892) should be around 1600 MeV. This is why we assign in Table V our 
1620 MeV state to the if (1630). It would be interesting for future experiments to check the 
quantum numbers of this state. With respect to the if*(1410) its assignment to the 2 3 5i 
state is not only excluded by its mass, but also by its decay modes. A possible interpretation 
of the low mass of this state due to the mixing with hybrids has been suggested in Ref. [74] 
and therefore seems to be excluded as a pure qq pair. 

The 2 + states are also in reasonable agreement with the theoretical predictions. This 
can be attributed to the understanding of the (IP) 2 ++ and (IP) 2 ++ isoscalar light-quark 
sector. The only difference is the mass of the strange quark in the interacting potential what 
makes clear the ns structure of if*(1430) and P|(1980). 

Strange mesons are the lighter ones where a mixing of the 1 + states with S = and 
S = 1 can occur (also the 2~ states) due to the fact that C parity is not a good quantum 
number. Although there is not theoretical consensus about the origin of such a mixing, 
in our model it is induced by the spin-orbit contributions of the OGE and confinement 
potentials proportional to S-. As a consequence the physical states would be given by 

| if*) = cos# I 1 + (S = 0)) - sin 6 | 1 + (S = 1)) 

| P x ) = sin 6> | 1+(S = 0)) + cos# | 1 + (S = 1)) (20) 

In the literature this mixing angle has been estimated by different methods. In Ref. [88] 
it was calculated from the ratio B(t — > z/Pi(1270))/P(r — > z/Pi(1400)), obtaining a value 
9 ~ +62°. In Ref. [74] an angle 9 ~ +45° was calculated from the pattern of the decay 
branching fractions of the two experimental states, if i (1270) and if i (1400). HQET [89] 
gives two possible mixing angles, 9 ~ +35.3° and 9 ~ —54.7°. We obtain a mixing angle of 
9 ~ 55.7° and therefore a mass for the physical states of 1352 MeV and 1414 MeV, improving 
the agreement with the experiment with respect to the unmixed masses. The mixing angle 
is close to the preferred value by the experimental decays and far from the results of Ref. 
[90] in a relativized quark model based only on OGE and confinement potentials. Their 
mixing angle 9 ~ +5° seems to be definitively excluded by the experimental decays. 

The same mixing occurs for the (2P) 1 + states, giving two resonances around 1.85 GeV, 
one would correspond to the Pi (1650) and the other to an unobserved state. Although 
this identification could seem inconsistent, the PDG states that the Pi(1650) entry contains 
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various peaks reported in partial-wave analysis in the 1600—1900 mass region. The election 
of a mass of 1650 MeV is done based on three measurements, one of 1650 ± 50 MeV and 
two around 1840 MeV. The same identification for 1^(1650) together with the prediction 
of a 2P tensor state 7^2(1850) has been done in Ref. [74] from the analysis of experimental 
data on strong decays. 

We have included in Table V a state given in the PDG as PJ 2 (1580) that has no equivalent 
in the qq spectrum. This state is clearly uncertain, it was reported in only one experimental 
work more than twenty years ago and has never been measured again. 

I. D and D s mesons 

Recently, a clear evidence for the existence of new open-charm mesonic states has been 
reported by three different collaborations, two states with charm-strange quark content and 
other with charm-light quark content. 

BaBar Collaboration reported a narrow state at 2316. 8±0. 5 MeV [62] called P*j(2317) . 
CLEO Collaboration [64,66,91] provided confirmation of the existence of this state and 
furthermore reported the observation of a new state called P s j(2463). In return, BaBar 
experiment also confirmed the existence of this state [63]. Finally both discoveries have been 
confirmed by the Belle Collaboration, not only in the analysis of inclusive e + e _ annihilation 
[67], but also in the exclusive B meson decays [68]. All the experimental observations are 
consistent with the assignment of P-wave states with spin-parity J p = + for the P*j(2317) 
and J p = 1 + for the P s j (2463). Belle Collaboration [59] has also reported the existence of 
a J p = + charm-light state with a mass of 2308±36 MeV. 

There are two intriguing aspects of the new D mesons. First of all they have a mass 
significantly smaller than the predictions of most QCD inspired quark potential models 
regarding these mesons as P-wave bound states, cs and cn respectively. Secondly they 
do not show the decay patterns favored by theoretical expectations. This has triggered a 
variety of articles either supporting the qq interpretation or presenting alternative hypothesis 
of exotic states. 

Our results are shown in Tables VI and VII. In Table VI we have only included those 
states reported by the PDG, while in Table VII we include all IP excited states. As can 
be seen the low-lying 0~ and l - D and D s states are perfectly reproduced. We obtain a 
theoretical state with the mass of the P*(2640) and spin-parity J p = 0~, and also a state 
with a similar mass to the P s j(2573) and J p = 2 + . However there is a discrepancy both 
with the J p = 1 + states and specially with the J p = + states. In the case of the 1 + states 
the same mixing as discussed in Eq. (20) between the 1 P\ and 3 Pi partial waves appears. 
As a consequence, for the charm-light sector we obtain two 1 + states with masses of 2454 
MeV and 2535 MeV, with a mixing angle of 43.5°, and for the charm-strange case 2543 MeV 
and 2571 MeV, with a mixing angle of 58.4°. Therefore, while we would find a candidate for 
the P s i(2536) and the Pi (2420), the recently measured P s j(2463) and its equivalent state 
in the charm-light sector (see Table VII) do not fit into the predictions of the nonrelativistic 
qq models [92]. 

The situation with the + states being in principle more discouraging is, however, similar 
to the problem observed for the light-scalar mesons (see Sec. Ill L), they hardly fit in a 
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qq scheme. As seen in Table VII the quark model predictions for the + D and D s are 
around 150 MeV above the experimental results. In a pure qq scheme the + states are 
influenced by the noncentral terms of the interaction, in particular they strongly depend 
on the spin-orbit force and therefore on the scalar/vector rate in the potential, controlled 
by a s . As this relation was fixed in the light-meson sector, where the potential has a much 
more involved structure (Goldstone boson exchanges), we wonder if the disagreement could 
be solved by modifying the scalar/vector rate of confinement. In Table XIII we present 
the results obtained by fixing a s to reproduce the experimental mass of one of the new 
measured D s mesons, the D* sJ (2317). Apart from obtaining a completely unusual value for 
the scalar/vector rate, a s = 0.46 [35,36], in doing this one observes how the situation of the 
charm-light sector does not improve, what makes evident the rather complicated situation 
appeared with the new measurements. 

Let us finally note that in a qq scheme the + states are obtained through an orbital 
angular momentum excitation, what could explain the large masses predicted. The + 
quantum numbers may also be obtained in the absence of orbital excitation considering 
more complex structures [93] which can help to understand the experimental situation. 

In the literature the new (0 + , 1 + ) states have triggered other alternative explanations 
as for example that they could be the members of the J p = (0 + , 1 + ) doublet predicted by 
the heavy quark effective theory (HQET) [94] or that these states together with the lowest 
(0~, 1~) constitute the chiral doublet of the (0+, 1+) HQET spin multiplet [95]. 



J. Charmonium and the new states 

We present in Table VIII our results for charmonium. Confinement parameters were 
fitted to describe the energy difference between the 5 = 1 ground state, J/ij)(lS), and 
its first radial excitation, ip(2S). One also observes a pretty good description of the first 
negative parity orbital excitation, -0(3770). The splitting among the Xcj states is correctly 
given both in order and magnitude. With respect to the S = sector, as happened for 
the light-quark systems, the orbital excitation is expected to be 600 MeV above the ground 
state, in complete agreement with the mass of the h c (lP). This state, although not yet 
confirmed is a clear candidate for being the 1 H resonance. There is also a set of L =even 
S = 1 states with an obvious correspondence to qq states although in some cases the S or 
D— wave identification is not unique, as seen on Table VIII for the 0(4415). 

Recently the Belle Collaboration has reported the observation of a narrow peak that has 
been interpreted as the 2S singlet charmonium state, the rj c (2S), 



M[ Vc (2S)} 



3654 ± 10 MeV Ref. [96] 
3622 ± 8 MeV Ref. [68] 
3630 ± 8 MeV Ref. [67] 



These masses are larger than the experimental value quoted by the 2002 PDG: M[r) c (2S)] = 
3594 ± 5 MeV and it was pointed out that cannot be easily explained in the framework of 
constituent quark models. The reason for that stems on the 2S hyperfine splitting (HFS) 
that would be smaller than the predicted for the IS* ones. In fact the theoretical predicted 
2S/1S HFS ratio is 
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R = 



AM 2 H S FS 
AMfi FS 



' 0.84 Ref. [97] 
< 0.67 Ref. [98] , 
k 0.60 Ref. [99] 



larger than the experimental value, R=0.273 if M[ri c (2S)] = 3654 MeV, R=0.547 for 3622 
MeV, and R=0.479 for 3630 MeV. Our result is M[rj c (2S)] = 3627 MeV, within the error 
bar of the last two Belle measurements, the ones obtained with higher statistics. Moreover 
the ratio 2S/1S HFS is found to be 0.537, in agreement with these last two experimental 
data. The reason for this agreement can be found in the radial structure of the confining 
potential that also influences the HFS, the linear confinement being not enough flexible to 
accommodate both excitations [46]. A similar agreement has also been obtained by other 
theoretical models including the effect of the open thresholds [47] (note that the result of 
Ref. [2] was far from the available experimental data at that time). 

Finally HFS is closely connected with the leptonic decay widths V (vector meson) — > 
e + e~. Although their absolute values depend on radiative and relativistic corrections, the 
ratios are a test of the wave functions at the origin, and still closely connected to HFS. We 
obtain T e + e - [ip{2S)]/Y e + e - [J/ip{lS)} =0.44, in good agreement with the experimental value 
0.41±0.07, what gives us confidence about the correct description of the rj c (2S). 

The most recently discovered charmonium state is the X(3872), reported by Belle [100] 
in the J/^n + 7i~ invariant mass distribution of the B ± — > K ± J/ip7i + 7i~ reaction, and con- 
firmed by the CDF collaboration at Fermilab [101]. Both experiments report a similar mass, 
3872.0 ± 0.6 ± 0.5 MeV in Ref. [100] and 3871.4 ± 0.7 ± 0.4 MeV in Ref. [101], very close to 
the D°D*° threshold (3871.5 ± 0.5 MeV). 

The proposed interpretations for the X(3872) include 1 3 L> 3 , 1 3 D 2 , 1 1 D 2 , 2 3 P 1 and 2 1 P 1 
cc states. None of them comfortably fit the observed properties of this state and therefore 
a considerable experimental uncertainty still remains. The results of our model for the 
possible quantum numbers are shown in Table IX. Although the well-established l 3 Di 
is well reproduced by our model the predicted ID candidates lie 70—80 MeV below its 
experimental mass whereas the 2P states are 40 MeV above. These results clearly show 
that the theoretical splitting in the 3 Dj multiplet is much smaller [~25 MeV] than the one 
necessary to correctly described the X(3872) and the \l/(3770) [~100 MeV]. Similar results 
have been found in other potential models which suggests that the X(3872) could present a 
more involved structure [102,103]. 



A pretty good description of the experimental states shown in Tables X and XI is ob- 
tained. Some caution is necessary with respect of some of the experimental data reported 
in Table XL The most important discrepancy observed between the data and our results 
arise from the rjb(lS), all the other states being described with similar accuracy to any other 
spectroscopic model designed to study a particular sector [46] . However, the result for this 
state is based in only one experimental work where only one event has been observed [45], 
and therefore this data needs confirmation. The uncertainty on this data can be easily un- 
derstood from the surprisingly large spin splitting that it would produce: my — m rjb = 160 



K. Bottomonium and open beauty states 



17 



MeV. Such an energy difference clearly spoils the evolution of the spin splittings with the 
mass of the quarks: m p — m n = 630 MeV, uik* — itir = 397 MeV, mo* — mo = 146 
MeV, m J/$ - m Vc = 117 MeV, m B * - m B = 46 MeV, m T - m Vb = 160 MeV. The S or D 
assignment of 1 high energy bb excitations is not conclusive. Let us mention again that 
with the noncentral terms of our interaction a correct description of the hyperfine splittings 
both in order and magnitude is obtained. For the open-beauty mesons (Table X) we find 
theoretical states corresponding to _B}(5732) and _B*j(5850), in both cases being J p = 2 + . 

L. J PC = 0++ states 

It is still not clear which are the members of the ++ nonet corresponding to the L = 
S — 1 nn and ss multiplets. There are too many ++ mesons observed in the region below 
2 GeV to be explained as qq states. There have been reported in the PDG two isovectors: 
a (980) and a (1450); five isoscalars: / (600), f (980), / ( 1370), / ( 1500) and / (1710); and 
three I = 1/2: 1^(1430), 1^(1950) and recently k(800). The quark model predicts the 
existence of one isovector, two isoscalars and two 1=1/2 states for each nonet. Our results 
are shown in Table XII. Using this table one can try to assign the physical states to ++ 
nonet members. 

Let us discuss each state separately. With respect to the isovector states, there is a 
candidate for the ao(980), the 3 Po member of the lowest 3 Pj isovector multiplet. The other 
candidate, the ao(1450), is predicted to be the scalar member of the 2 3 Pj excited isovector 
multiplet. This reinforces the predictions of the quark model, the spin-orbit force making 
lighter the J = states with respect to the J = 2. The assignment of the ao(1450) as 
the scalar member of the lowest 3 Pj multiplet [104] would contradict this idea, because the 
a2(1320) is well established as a qq pair. The same behavior is evident in the cc and bb 
spectra, making impossible to describe the ao(1450) as a member of the lowest 3 Pj isovector 
multiplet without spoiling the description of heavy-quark multiplets. However, in spite of 
the correct description of the mass of the a (980), the model predicts a pure light-quark 
content, what seems to contradict some experimental observations [105]. The a (1450) is 
predicted to be also a pure light-quark structure obtaining a mass somewhat higher than 
the experiment. 

In the case of the isoscalar states, we find a candidate for the / (600) with a strangeness 
content around 10%. There are no / = states with a mass close to 1 or 1.5 GeV, which 
would correspond to the /o(980) and the /o(1500), and they cannot be found for any com- 
bination of the model parameters. It seems that a different structure rather than a naive qq 
pair is needed. In particular, the /o(1500) is a clear candidate for the lightest glueball [106], 
while the / (980) has been suggested as a possible four-quark state [107,108] what would 
make it compatible with the similar branching ratios observed for the J/ip — > /o(98O)0 and 
J/ijj — > /o(980)u; decays [105]. Concerning the /o(1370) (which may actually correspond to 
two different states [109]) we obtain two states around this energy, the heavier one with a 
dominant nonstrange content which favors its assignment to the /o(1370); the other with 
a high ss content without having an experimental partner. Let us however remember that 
in this energy region there is a state, the / 2 (1430), that does not fit into the isoscalar 2 ++ 
sector and has been recently suggested as a possible ++ state [83]. Finally a dominant 
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nn state corresponding to the /o(1710) is obtained. Our results concluding that /o(1370), 
/o(1500) and /o(1710) are dominantly nn, non qq, and nn respectively differ from the con- 
clusion of Refs. [106,110] obtaining that / (1710) is dominantly ss and are also in contrast 
to the predictions of Ref. [Ill] which prefers to assign / (1500) to an ss state and / (1710) 
to a glueball. This makes clear the complicated situation in the scalar sector with several 
alternative interpretations of the observed states. The study of radiative transitions and 
two photon decay widths should help to understand the flavor mixing and the nature of the 
1 = scalar sector. We obtain two states around 1.9 GeV, a 2 3 P state with a dominant 
ss content and a 4 3 P with a dominant nn content. Being the / (2020) an experimentally 
known nn meson, its identification with the 4 3 P state is clear. Therefore, we find an un- 
observed ss scalar with a mass around 1.9 GeV as has also been suggested in Ref. [74]. 
Finally, although for consistency not given in the tables, we find a candidate for the 3 3 Po 
state / (2200), experimentally identified as an ss state [112], with an energy of 2212 MeV. 

Concerning the 1=1/2 sector, as a consequence of the larger mass of the strange quark 
as compared to the light ones, our model always predicts a mass for the lowest ++ state 
200 MeV greater than the ao(980) mass. Therefore, being the ao(980) the member of the 
lowest isovector scalar multiplet, the k(800) cannot be explained as a qq pair. We find a 
candidate for the -£To(1430) although with a smaller mass. 

Our results indicate that the light-scalar sector cannot be described in a pure qq scheme 
and more complicated structures or mixing with multiquark system seems to be needed. 
Concerning the /o(980) and the k(800) our conclusions are very similar to those obtained 
in Ref. [113] using the extended Nambu-Jona-Lasinio model in an improved ladder approx- 
imation of the Bethe-Salpeter equation. This seems to indicate that relativistic corrections 
would not improve the situation and the conclusions remain model independent. 

One finds in the literature several alternative explanations to understand the rather 
complicated scenario of the scalar mesons. An earlier attempt to link the understanding 
of the NN interactions with meson spectroscopy was done based on the Julich potential 
model [114]. The structure of the scalar mesons ao(980) and /o(980) was investigated in 
the framework of a meson exchange model for ttti and nn scattering. The KK interaction 
generated by the vector-meson exchange, which for isospin / = is strong enough to generate 
a bound state is much weaker for / = 1, making a degeneracy of a (980) and / (980) 
impossible, as found in our model. Although both scalar mesons result from the coupling 
to the KK channel explaining in a natural way their similar properties, the underlying 
structure obtained was, however, quite different. Whereas the / (980) appears to be a KK 
bound state the a (980) was found to be a dynamically generated threshold effect. 

In a different fashion within the quark model the same problem was illustrated in Ref. 
[115]. The bare mass used for the nn pair is much larger than the ao(980) and /o(980) 
experimental masses. It is the effect of the two-pseudoscalar meson thresholds the responsible 
for the substantial shift to a lower mass than what is naively expected from the qq component 
alone. This gives rise to an important KK and nn' components in the a (980) and KK, nn, 
n'n' and nn' in the / (980). In particular for the a (980) they obtain the KK component 
to be dominant near the peak, being about 4 — 5 times larger than the qq component. 
A similar conclusion, that the description of the ao(980) and /o(980) requires from more 
complex structures, is also obtained from our analysis. The absence of the /o(1500) in our 
qq scheme make also contact with the indication that this state could correspond to the 
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lightest scalar glueball [116]. 

A similar problem as the one observed with the light scalars appeared in the open charm 
sector. As we have already discussed, the two recently measured + states do not fit into 
a qq description. Possible alternatives to understand their masses, as for example being 
DK molecules [117] or tetraquarks [107], have been suggested. Using the same interacting 
hamiltonian presented in this work an estimation of the lowest scalar open charm tetraquark 
has been done in Ref . [108] . The variational estimation was performed under the assumption 
that internal orbital angular momentum does not give an important contribution, what has 
been strictly tested in the case of heavy-light tetraquarks [118], but it is well known to 
influence light-quark systems. Having in mind this precaution a mass of 2389 MeV was 
obtained for the light open-charm [(ns)(nc)] tetraquark that could very well correspond to 
the _D*j (23 17). A consistent calculation of four light-quark structures seems to be advocated 
for a full understanding of the scalar sector. 

IV. SUMMARY 

We have performed an exhaustive study of the meson spectra from the light nn states to 
the bb mesons within the same model. The quark-quark interaction takes into account QCD 
perturbative effects by means of the one-gluon-exchange potential and the most important 
nonperturbative effects through the hypothesis of a screened confinement and the sponta- 
neous breaking of chiral symmetry. The model incorporates in a natural way the nn <-> ss 
isoscalar mixing due to the Vk potential, and the 1 + S = and S = 1 mixing caused by 
the spin-orbit force. An arbitrary rate of scalar/vector confinement has been used, finding 
evidence of a strong scalar component. Annihilation effects have not been taken into ac- 
count although their contribution seems to be important for the description of the l +_ light 
isoscalar ground state. We have obtained a reasonable description of most part of the well 
established qq states for all flavors. The success of the model allowed us to make predictions 
with respect to those states whose quantum numbers, existence or nature is under debate. 

In the light-isoscalar sector our results support the speculation pointed out by the PDG 
about the possible nonexistence of the /i(1420). We do not find a theoretical 2 ++ state 
corresponding to the /2(1430), in agreement with recent experiments that opened the pos- 
sibility of a different assignment of quantum numbers for this meson being a ++ state. 
Our model predicts a scalar meson in this energy region without an experimental partner. 
Besides, we have only found one qq state in the 1.6 GeV region, which seems to indicate 
that either the /2(1565) and the /2(1640) are the same state as suggested by the PDG, or 
that the /2(1565) goes beyond the naive qq structure. In the light-scalar sector it seems very 
difficult to accommodate the / (980), the k(800) and the / (1500) in a qq scheme. In the 
light-strange sector we do not find a qq state to be identified with the if* (1410), favoring 
its possible hybrid structure. The same situation occurs for the K 2 (1580), although in this 
case the poor experimental data do not assure its existence. 

Concerning the flavor content of the 7/(1440) our model indicates that it is a dominant 
nn state, with a probability of 70.3 %. The 03(1875) is compatible with an ss content but 
there should exist an unobserved nn state in the same energy region. We have also found 
evidence of the existence of a l + ~ light-isoscalar meson with a dominant ss content and a 
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mass around 1.97 GeV, and an isoscalar ss 2 + ~ state at 1.85 GeV. Finally, our model assigns 
the quantum numbers 1~ to A" (1630). 

For the heavy-quark sector the experimental situation is changing very fast. New ex- 
periments and reanalysis of old data are being done and new states being discovered. Some 
of them fit nicely in a qq scheme, but others are impossible to accommodate. In the open 
charm sector, the D*(2640) and D s j(2573) are compatible with cn and cs mesons with 
quantum numbers 0~ and 2 + , respectively. However, the recently discovered + states, the 
-D*j(2317) and _Dj(2308), seem to have a completely different structure. We have argued a 
possible explanation to describe their low masses based in a tetraquark structure that may 
have positive parity without orbital excitation. Finally, there is an obvious identification of 
a qq state for the h c (lP) with quantum numbers 1 H . 

We consider that this type of study based on models whose parameters are constrained in 
the description of other low-energy systems should be a complementary useful tool to deepen 
the understanding of the meson spectra. The next step in this effort to a comprehensive 
description of the new data concerning the meson spectra should be the analysis of the 
electroweak and strong decays of mesons that will be the subject of a future publication 
[119]. 
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TABLES 



TABLE I. Quark model parameters. 
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TABLE II. Energies (in MeV) obtained when solving the Schrodinger equation for two different 
partial waves when the different noncentral terms are switched on. 



Potential 

+V c SO if ij ) 



(1=1) 
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1598 
1474 
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TABLE III. Masses, in MeV, of / = light-quark mesons up to 2 GeV. QM denotes the results 
of the present model and Flavor stands for the dominant component of the flavor wave function. 
Experimental data, PDG, are taken from Ref. [45]. In the second column we denote by a question 
mark those states whose existence is not clear. In the third column, if there are several candidates 
for an experimental state, we underline our preferred assignment. See text for details. 
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TABLE IV. Masses, in MeV, of I = 1 light-quark mesons up to 2 GeV. QM denotes the 
results of the present model. Experimental data (PDG) are taken from Ref. [45]. We denote by a 
dagger those PDG states whose masses are not explicitly given. We use for them the more recent 
experimental data. 
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TABLE V. Masses, in MeV, of the light-strange mesons up to 2 GeV. QM denotes the results of 
the present model. Experimental data (PDG) are taken from Ref. [45]. We denote by a dagger those 
PDG states whose masses are not explicitly given. We take for them the more recent experimental 
data. The spin is indicated because C parity is not well defined. In those cases where the PDG 
does not give the J p quantum numbers and we find a candidate for the state, we quote between 
square brackets our predictions. QM(mixed) denotes the mass of the states after being mixed 
according to Eq. (20). 
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1650±50 

1425.6±1.5 
1973±26 
« 1580 f 

1773±8 
1816±13 

1776±7 
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TABLE VI. Same as Table V for D and D s mesons. 



Meson 


{nL) J p 


Spin 


State 


QM 


QM(mixed) 


PDG 




(15)0- 





D 


1883 




1867.7±0.5 




(15) 1- 


1 


D* 


2010 




2008.9±0.5 


D 


(1P) 1+ 





P>i(2420) 


2492 


2454 


2425±4 




(IP) 2+ 


1 


D|(2460) 


2502 




2459±4 




(25 s ) ? ? ro — l 

\AO ) . [U J 


o 


n* (2640 N l 


2642 




2637±7 




(15)0- 





D s 


1981 




1968.5±0.6 




(15) 1- 


1 


D* s 


2112 




2112.4±0.7 


D s 


(1P) 0+ 


1 


D* sJ (2317) 


2469 




2317.4±0.9 




(IP) 1+ 





^(2460) 


2550 


2543 


2459.3±1.3 




(IP) 1+ 


1 


P» si (2536) 


2563 


2571 


2535.3±0.6 




(IP) ? ? [2+] 


1 


D sJ {2573) 


2585 




2572.4±1.5 



TABLE VII. Masses, in MeV, of the first positive parity D and D s mesons compared to recently 
measured experimental data. QM denotes the results of the present model. QM(mixed) indicates 
the mass of the states after being mixed according to Eq. (20). 



Meson 


(nL) J p 


(IP) 0+ 


(IP) 1+ 


(IP) 1+ 


(IP) 2+ 


D 


PDG [45] 
FOCUS [57,58] 
Belle [59] 
CLEO [60] 
DELPHI [61] 


~ 2420 
2308±36 


2427±36 
2461±51 
2470±58 


2425±4 
2421±2 


2459±4 
2468±2 
2461±4 




QM 
QM (mixed) 


2436 


2496 
2535 


2492 
2454 


2502 


D s 


PDG [45] 
FOCUS [57,58] 
BaBar [62,63] 
CLEO [64-66] 
Belle [67,68] 


2317.4±0.9 

2317.3±0.9 
2318. 1±1. 2 
2319.8±2.1 


2459.3±1.3 

2458.0±1.4 
2463.1±2.0 
2456.5±1.8 


2535.3±0.6 
2535.1±0.3 


2572.4±1.5 
2567.3±1.4 




QM 
QM(mixed) 


2469 


2563 
2571 


2550 
2543 


2585 
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TABLE VIII. Same as Table III for charmonium. 



(nL) J PC 


State 


QM 


PDG 


(15) o-+ 


Vc(lS) 


2990 


2979.6±1.2 


(15) 1- 


J/V(15) 


3097 


3096.916±0.011 


(1P) 0++ 


Xco(lP) 


3436 


3415.19±0.34 


(IP) 1++ 


Xcl(lP) 


3494 


3510.59±0.10 


(IP) 2++ 


Xc2(lP) 


3526 


3556.26±0.11 


P) ? ?? [1+-] 


/i c (lP) 


3507 


3526.21±0.25 


(25)0-+ 


r? c (25) 


3627 


3654±10 


(25)1- 


^(25) 


3685 


3686.093±0.034 


(LD)l- 


V>(3770) 


3775 


3770.0±2.4 


(1£>)2- 


V>(3836) 


3790 


3836±13 


(35)1- 


V>(4040) 


4050 


4040±10 


(2L>) 1- 


V>(4160) 


4103 


4159±20 


(45)1- 
(3P>) 1 


V>(4415) 


" 4307 " 
4341 


4415±6 



TABLE IX. Masses predicted by our model for the qq states compatible with the X(3872) 



1 3 D 3 


1 3 D 2 




l l D 2 


2 3 P 1 


2 x Pi 


3802 


3790 




3793 


3913 


3924 






TABLE X. 


Same as 


Table V for B, B s and B c 


mesons. 








(nL) J p 


Spin 


State 


QM 


PDG 


B 




(15)0- 
(15)1- 
(1P) ? ? [2+] 



1 
1 


B 

B* 
P}(5732) 


5281 
5321 
5790 


5279.2±0.5 
5325.0±0.6 
5698±8 


B s 




(15) 0" 
(15)1- 
(1P) ? ? [2+] 



1 
1 


B s 

B* s 
B* sJ (5850) 


5355 
5400 
5855 


5369.6±2.4 
5416.6±3.5 
5853±15 


B c 




(15)0- 





B c 


6277 


6400±410 



31 



TABLE XI. Same as Table III for bottomonium. We denote by an asterisk a experimental 
state recently reported in Ref. [69]. 



(nL) J PC 


State 


QM 


PDG 


(15)0-+ 


Vb(lS) 


9454 


9300±28 


(15)1- 


T(15) 


9505 


9460.30±0.26 


(1P) 0++ 


Xbo(lP) 


9855 


9859.9±1.0 


(IP) 1++ 


Xbi(lP) 


9875 


9892.7±0.6 


(IP) 2++ 


Xb2(lP) 


9887 


9912.6±0.5 


(25)1- 


T(25) 


10013 


10023.26±0.31 


(10)2- 


T(LD 2 )* 


10119 


10162.2±1.6* 


(2P) 0++ 


Xfeo(2P) 


10212 


10232. 1±0.6 


(2P) 1++ 


X6i(2P) 


10227 


10255. 2±0. 5 


(2P) 2++ 


X6 2 (2P) 


10237 


10268. 5±0.4 


(35)1- 


T(35) 


10335 


10355.2±3.5 


(45)1" 


T(45) 


10577 


10580.0±3.5 


(55)1-- 
(AD) 1— 


T(10860) 


' 10770 " 
10803 


10865±8 


(65) 1- 
(5D) 1 


T(11020) 


' 10927 " 
10953 


11019±8 



TABLE XII. Same as Table III for the light-scalar mesons. We have included the rc(800), whose 
isospin is not well determined being preferred 1 = 1/2 [70]. 





(nL) J PC 


State 


QM 


Flavor 


PDG 


I = 1 


(IP) 0++ 
(2P) 0++ 


a (980) 
a (1450) 


984 
1587 


(nn) 
(nn) 


984.7±1.2 
1474±19 


7 = 




(IP) 0++ 
R0++ 
(2P) 0++ 
(IP) [0++] 

(-)o++ 

(3P) 0++ 
' (4P) 0++ " 
(2P) 0++ 


/o(600) 
/o(980) 
/o(1370) 

/o(1500) 
/o(1710) 

/ (2020) 


413 

1395 
1340 

1754 
' 1880 " 
1894 




(nn) 

(nn) 
(ss) 

(nn) 
(nn) 
(ss) 


400-1200 

980±10 
1200-1500 

1507±5 
1714±5 

1992±16 


I = 1/2 


(-)0+ 
(IP) 0+ 
(2P) 0+ 


k(800) 
K^(1430) 
K^(1950) 


1213 
1768 


(ns) 
(ns) 


w 800 
1412±6 
1945±22 
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TABLE XIII. Masses, in MeV, of the first positive parity D and D s mesons with a s = 0.46. 
They are already mixed according to Eq. (20). 



Meson 


nL) J p 


State 


QM (mixed) 




(IP) 0+ 


££,(2317) 


2317 


n 

Us 


(IP) 1 + 




2482 




(IP) 1+ 


P s i(2536) 


2574 




(IP) 2+ 


P St7 (2573) 


2633 




(IP) 0+ 


2308±36 


2134 


D 


(IP) 1+ 


Pi (2420) 


2354 




(IP) 1+ 


2427±36 


2524 




(IP) 2+ 


P|(2460) 


2588 
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FIGURES 



FIG. 1. Effective scale-dependent strong coupling constant a s given in Eq. (14) as a function 
of momentum. We plot by the solid line our parametrization. Dots and triangles are the exper- 
imental results of Refs. [24] and [25], respectively. For comparison we plot by a dashed line the 
parametrization obtained in Ref. [18] using A = 0.2 GeV. 

FIG. 2. Masses of the aj triplet members as a function of the scalar/vector rate confinement, 
a s . The solid line denotes the ao mass, the dashed that of a\ and the dashed-dotted stands for the 
mass of the a-z meson. A vertical solid line indicates the value chosen for a s . 
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